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The one-dimensional penetrable-square-well fluid is studied using both analytical tools and spe- 
cialized Monte Carlo simulations. The model consists of a penetrable core characterized by a finite 
repulsive energy combined with a short-range attractive well. This is a many-body one-dimensional 
problem, lacking an exact analytical solution, for which the usual van Hove theorem on the absence 
of phase transition does not apply. We determine a high-penetrability approximation complementing 
a similar low-penetrability approximation presented in previous work. This is shown to be equivalent 
to the usual Debye-Hiickel theory for simple charged fluids for which the virial and energy routes 
are identical. The internal thermodynamic consistency with the compressibility route and the va- 
lidity of the approximation in describing the radial distribution function is assessed by a comparison 
against numerical simulations. The Fisher- Widom line separating the oscillatory and monotonic 
large-distance behavior of the radial distribution function is computed within the high-penetrability 
approximation and compared with the opposite regime, thus providing a strong indication of the 
location of the line in all possible regimes. The high-penetrability approximation predicts the ex- 
istence of a critical point and a spinodal line, but this occurs outside the applicability domain of 
the theory. We investigate the possibility of a fluid-fluid transition by Gibbs ensemble Monte Carlo 
techniques, not finding any evidence of such a transition. Additional analytical arguments are given 
to support this claim. Finally, we find a clustering transition when Ruelle's stability criterion is 
not fulfilled. The consequences of these findings on the three-dimensional phase diagrams are also 
discussed. 



I. INTRODUCTION 

Recent advances in chemical synthesis have unveiled 
more and more the importance of soft-matter systems, 
such as dispersions of colloidal particles, polymers, and 
their combinations. Besides their practical interest, these 
new developments have opened up new theoretical av- 
enues in (at least) two instances. Firstly, it is possi- 
ble to experimentally fine-tune the details of interactions 
(range, strength, . . . ) , making these systems a unique 
laboratory for testing highly simplified models within an 
effective-interaction approach where the microscopic de- 
grees of freedom are integrated out in favor of renormal- 
ized macroparticle interactions. Secondly, they offer the 
possibility of exploring new types of equilibrium phase 
behaviors not present in the simple-fluid paradigm. 

As early as in 1989, Marquest and Witten^ suggested 
that the experimentally observed crystallization in some 
copolymer micellar systems could be rationalized on the 
basis of a bounded interaction, that is, an interaction 
that does not diverge at the origin. Successive theoreti- 



cal work showed that this class of bounded or ultra-soft 
potentials naturally arises as effective interactions be- 
tween the centers of mass of many soft and flexible macro- 
molecules, such as polymer chains, dendrimers, star poly- 
mers, etc. (see, e.g., Ref. 2 for a reference on the sub- 
ject). Two well-studied cases belonging to the above 
class are the Gaussian core model (GCM) introduced by 
Stillinger'^ and the penetrable-sphere (PS) model intro- 
duced in Refs. 1,4, whose freezing transition turns out to 
display rather exotic features with no analogue in the 
atomistic fluid realm. 

In the present paper, we shall consider a close relative 
of the PS model, first introduced in Ref. 5, denoted as 
the penetrable-square-well (PSW) model, where a short- 
range attractive tail is added to the PS model just outside 
the core region. In the limit of infinite repulsive energy, 
the PS and PSW models reduce to the usual hard-sphere 
(HS) and square-well (SW) models, respectively. 

An additional interesting feature common to both PS 
and PSW, as well as to all bounded potentials, is the 
fact that, even when confined to one-dimensional sys- 



tcms, they may exhibit a non-trivial phase diagram due 
to the penetrability which prevents an exact analytical 
solution. 

This is because particles cannot be lined up on a line 
with a well-defined and fixed ordering in view of the pos- 
sibility of reciprocal interpenetration (with some posi- 
tive energy cost), thus lacking an essential ingredient 
allowing for an exact solution in the respective hard- 
core counterparts (HS and SW). It is then particularly 
useful to discuss some motivated approximations in the 
one-dimensional model which can then be benchmarked 
against numerical simulations and subsequently exploited 
in the much more complex three-dimensional case. 

The aim of the present paper is to complete a study 
on the one-dimensional PSW model started in Refs. 5,6, 
as well as the general results presented in Ref. 5 which 
are particularly relevant in the present context. In the 
first paper of the series,^ we introduced the model and 
discussed the range of stability in terms of the attractive 
versus repulsive energy scale. We provided, in addition, 
exact analytical results in the low-density limit (second 
order in the radial distribution function and fourth or- 
der in the virial expansion) and a detailed study of the 
Percus-Yevick and hypernetted-chain integral equations. 
These were used in the following papcr^ to propose a 
low-penetrability approximation at finite density which 
was then tested against numerical simulation. This low- 
penetrability approximation is expected to break down 
in the opposite regime, namely when temperatures and 
densities are such that particles easily interpenetrate each 
other. In the present paper, we address this latter regime 
by proposing a complementary approximation (the high- 
penetrability approximation) and discussing its range of 
validity and the relationship with the low-penetrability 
regime. Note that a similar matching of the low- and 
high-penetrability approximation has already been car- 
ried out by two of the present authors in the framework of 
the PS model. ^'® It turns out that the high-penetrability 
approximation in the context of bounded potentials coin- 
cides with the linearized Debye-Hiickel classical approx- 
imation originally introduced in the framework of the 
Coulomb potential.^ It has been recently shown^° that 
two of the three standard routes to thermodynamics (the 
energy and the virial routes) are automatically consis- 
tent within the linearized Debye-Hiickel approximation, 
for any potential and dimensionality. This means that 
a deviation from the third standard route to thermody- 
namics (the compressibility route) can be exploited to 
assess the degree of reliability of the high-penetrability 
(or linearized Debye-Hiickel) approximation. This is in- 
deed discussed in the present paper, where we also discuss 
the full hierarchy of approximations ranging from the 
full Debye-Hiickel approximation to the simplest mean- 
spherical approximation. 

In view of the boundness of the potential, the usual van 
Hove no-go theorem^ ^'^^ on the absence of phase tran- 
sitions in certain one-dimensional fluids does not hold. 
It is then natural to ask whether a phase transition oc- 



curs in the one-dimensional PSW fluid by noting that 
the addition of an attractive tail to the pair potential 
of the PS model extends the question to the fluid-fluid 
transition, in addition to the fluid-solid transition pos- 
sible even within the PS model. In the present paper 
we confine our attention to the fluid-fluid case only and 
discuss this possibility using both analytical arguments 
and state-of-the-art numerical simulations. ^'^"^^ Our re- 
sults are compatible with the absence of such a transi- 
tion, as we shall discuss. This is also supported by recent 
analytical resultS"'^'' using a methodology devised for one- 
dimensional models with long-range interactions.^® We 
discuss possible reasons for this and a plausible scenario 
for the three-dimensional case. 

Finally, we note that the approach to a critical point 
is frequently anticipated by the so-called Fisher- Widom 
line^^ marking the borderline between a region with oscil- 
latory behavior in the long-range domain of the correla- 
tion function (above the Fisher- Widom line) and a region 
of exponential decay. We discuss the location of this line 
within the high-penetrability approximation and again 
the matching of this result with that stemming from the 
low-penetrability approximation. 

The structure of the paper is as follows. We deflne 
the PSW model in Sec. II. We then construct the high- 
penetrability approximation in Sec. HI and in Sec. IV we 
discuss some approximations related to it. Section V con- 
tains a discussion on the routes to thermodynamics, as 
predicted by the high-penetrability approximation. The 
structure predicted by the approximation is compared 
with Monte Carlo data in Sec. VI. The Fisher-Widom 
line and the possibility of a fluid-fluid transition are dis- 
cussed in Sec. VII. The paper ends with some concluding 
remarks in Sec. VIII. 



II. THE PENETRABLE-SQUARE- WELL (PSW) 
MODEL 



The penetrable-square- well (PSW) model is defined by 
the following pair potential^'^ 



0PSW (r) 



+£,., r < a, 
-Ea, (7 <r < a - 
0, r>a + A, 



A, 



(2.1) 



where A is the well width and e^ and Eq arc two pos- 
itive constants accounting for the repulsive and attrac- 
tive parts of the potential, respectively. In the following, 
we shall restrict our analysis to the case A/a < 1 and 
er > 2eai where we know the one-dimensional model to 
be stable with a well defined thermodynamic limit. ^ It 
is shown in Appendix A that, more generally, the one- 
dimensional PSW model is guaranteed to be stable if 
Cr > 2{£ + l)ea, where £ is the integer part of A/a. For 
lower values of Cr the model may or may not be stable 
and we will come back to this point in Section VII. 
An important role in the following is played by the 



corresponding Mayer function 

/psw (r) = 7r/HS (r) + 7a [6 (r - cr) - 6 (r - cr - A)] , 

(2.2) 



where 7^ = 1 



-I3€^ 



is the parameter measuring the 



degree of penetrability varying between (free penetra- 
bihty) and 1 (impenetrability), while 7^ = e'^'^" — 1 > 
measures the strength of the well depth. Here /3 = l/fcsT 
with T the temperature and k^ the Boltzmann constant, 
/hs(^) — — ©(c — r) is the Mayer function for hard 
spheres, and Q{r) is the Heaviside step function. 

A detailed discussion of the limiting cases of the PSW 
model can be found in Ref. 5. Here we merely note that 
the PSW Mayer function /psw('') is immediately related 
to the usual SW Mayer function by 



fpsw{r) = 7r/sw (r) , 



(2.3) 



/sw (r) = /hs (r) + 7 [9 (r - ct) - 9 (r - cr - A)] , 

(2.4) 
where we have introduced the ratio 7 ~ "fa/jr- At a 
given value of er/ea, 7 increases quasi-linearly with e^*^", 
its minimum value 7 = ea/^r corresponding to (3ea ^ 0. 



III. A HIGH-PENETRABILITY 
APPROXIMATION (HPA) 

In Ref. 6 we discussed a low-penetrability approxima- 
tion (LPA) to the PSW model. Within this approxima- 
tion, one assumes 1— 7r = e~^'^'' ^ 1 so that the repulsive 
barrier e^ is sufficiently higher than the thermal energy 
ksT, penetrability is small, and the system is almost a 
hard-core one. The advantage of this theoretical scheme 
is that one can use the general recipe leading to the exact 
solution for the one-dimensional SW problem — in fact, 
valid for any potential with a hard core and short-range 
attractions — and perform some ad hoc adjustments to 
ensure that some basic physical conditions on the radial 
distribution function g{r) are satisfied. Comparison with 
Monte Carlo (MC) simulations showed a good behavior 
of the LPA even for f3er — 2 {1 ~ jr — 0.14), provided 
the density was moderate (pa < 0.5). 

The opposite limit 7^ ^ 1 is also inherently interest- 
ing for several reasons. From a physical viewpoint this 
amounts to starting from the ideal gas limit 7^ — >■ (one 
of the common reference systems for simple fluids) and 
progressively building up interactions by increasing 7^. 
An additional mathematical advantage stems from the 
simple observation^'^'^" that in the (exact) cluster expan- 
sion of g{r) only the dominant chain diagrams need to be 
retained at all orders, thus leading to the possibility of an 
exact summation of those leading contributions. As we 
shall see shortly, this is in fact a procedure known as the 
Debye-Hiickel approximation in the context of charged 
fiuids.9'10 

Our main goal is the computation of the cavity func- 
tion y{r) = e^'^^'^'g{r), from which one can immedi- 
ately compute the radial distribution function (RDF) 



g{r) = y{r)[l + f{r)]. In the PSW case one then has 
from Eq. (2.3) 

{(1 -lr)y{r), r < 1, 
(l+7a)2/(r), l<r<A, (3.1) 

y(r), r > X. 

where A = 1 + A and, in conformity with previous 
work,^'^ we have redefined all lengths in units of a so 
we set (T = 1 in most of the following equations. 

As shown in Ref. 20 for the PS case, the exact form 
of the PSW cavity function in the limit 71- — > at finite 
7 = 7a /7r and pj,. is 



y{r) = 1+ "frW (r) , 



(3.2) 



where the function w{r) is defined through its Fourier 
transform 



wik) 



Plr 



JSWV^) 



1 - Plrfsw{k)' 



(3.3) 



fsw{k) being the Fourier transform of /sw(^)- Note that 
in the limit 7^ ^- one has jr ~ f^Sr and 7 w Sa/^r- 

Generalizing an analogous approximation in the con- 
text of the PS model, ^'® our high-penetrability approxi- 
mation (HPA) consists of assuming Eqs. (3.2) and (3.3) 
for finite, but small, values of 7^. It is worth noting that 
the combination of the two expressions (3.2) and (3.3) de- 
fines what is usually referred to, in a different context, as 
the linearized Debye-Hiickel (LDH) approximation, ^'^° 
and this will be further elaborated below. 

Equations (3.2) and (3.3) hold for any dimensional- 
ity. In the specific one-dimensional case, and taking into 
account Eq. (2.4), we have 

/•CXD 

/sw(fc) = 2 / drcos{kr)fsw{r) 
Ja 
2 

= -T[(l + 7)sinfc-7sinA/i:] . (3.4) 
k 

The function 'w{r) can be numerically evaluated in real 
space by Fourier inversion as 



w{r) 



Plr 
■K 



dk cos(fcr) — 
1 



fswi'^) 



(3.5) 



P7r/sw(fc) 

An explicit expression for the density expansion of w{r) 
within the HPA is reported in Eqs. (B7) and (B8) of 
Appendix B, where the radius of convergence of the ex- 
pansion is also analyzed. 

From Eqs. (2.3), (3.1), and (3.2) the total correlation 
function, h{r) = g{r) — 1, within the HPA is easily ob- 
tained as 

h (r) = "frw (r) [1 + 7r/sw (r)] + 7r/sw {r) , (3.6) 
or in Fourier space, 

/sw (fc) 



h{k) = 7r 



+Plr 



1 ~ Plrfsw{k) 

°° dfc'/|w(fc')/sw(|fc-fc'|) 



27r 1 - pjrfsw {k') 



(3.7) 



From this equation it is straightforward to get the struc- 
ture factor 

S{k) = l + ph{k) 
1 



7r/sw (k) or, in real space, 



1 - Plrfsw (k) 

2 3 / dfc'/|w(fc')/sw(|fc-fc'|) 



+^^^/ ^ 



1 - PlrfsW {k') 



(3.8) 



and the Fourier transform of the direct correlation func- 
tion c(fc) = h{k)/S{k). The zero wavenumber value of 
the structure factor is 



S{0) 



1 



l + 2p7,(l-7A) 



-p'l'r r - 



fsw\k) 



h T^ l~ pjrfsW (fc) 

(3.9) 
where we have taken into account that /gw (0) = ^2(1 — 
7A). This completes the calculation of the correlation 
functions within the HPA. 



IV. APPROXIMATIONS RELATED TO THE 
HPA 

As anticipated in Section III, the HPA is the exact 
equivalent to the well-known LDH approximation, which 
is widely used in the context of charged simple fluids.^ 
The latter is actually an intermediate step of a hierarchy 
of successive approximations ranging from the simplest 
mean-spherical approximation (MSA) to the full non- 
linear version of the Debye-Hiickel approximation (see 
Ref. [11] in Ref. 10 for a discussion on this point). For 
the PSW model — and more generally for any bounded 
potential — the LDH approximation (the HPA in the 
present language) is particularly relevant in view of the 
fact that one can make /psw('') arbitrarily small by let- 
ting 7r — >■ 0, thus justifying the approximation of neglect- 
ing non-chain diagrams. It is then interesting to check 
the performance of the other approximations included in 
the aforementioned class, which will be translated in the 
present context for simplicity. 

On top of the hierarchy of approximations there is the 
non-linear HPA (nlHPA) 



yir) 



o7rtu(r) 



(4.1) 



which is equivalent to the non-linear Debye-Hiickel ap- 
proximation, as remarked. The HPA, Eq. (3.2), is ob- 
tained upon linearizing the exponential, an approxima- 
tion valid again in the limit 7^ <S 1. An additional 
approximation — denoted here as the modified HPA 
(mHPA) — can be considered with the help of Eq. (3.6) 
by neglecting the quadratic term in 7^. This yields 



^i (r) == 7r [w(r) + /sw (r)] 



(4.2) 



This is equivalent to keeping only the first term on 
the right-hand side of Eq. (3.7), which implies c(fc) — 



c{r) = fpswir) 



(4.3) 



The lowest rank in the hierarchy is occupied by the MSA, 
which is obtained from Eq. (4.3) upon linearization of the 
Mayer function /psw('')i 



c{r) 



-'I3(f>psw (r) 



(4.4) 



Since w{r) is a convolution, it must be continuous at 
r = 1 and r = X. It follows that the approximations 
with a continuous cavity function at r = 1 and r — X are 
nlHPA and HPA. For instance, in the mHPA (4.2) the 
cavity function is y{r) = 1-1- ^rw{r)/[l + jrfswi'r)], so 
that one has y(l-)-y(l+) = 72w(l)(l + 7)(l-7,,)-i(l + 
77r)-i andy(A+) -y{X-) = -f^w{X)j{l + j-fr^^ 

It has been shown in Ref. 10 that the virial and energy 
routes to thermodynamics (to be discussed below) are 
consistent one another within the HPA, for any poten- 
tial and any dimensionality. A similar statement holds 
true for soft potentials within the MSA.^^ This clearly 
includes the PSW potential in both cases. 

It is interesting to make contact with previous work 
carried out by Likos et al.^^'^'^ on a general class of un- 
bounded potentials which are free of attractive parts, 
thus resulting particularly useful in the context of the 
fiuid-solid transition. In Ref. 22, the MSA given in Eq. 
(4.4) along with the spinodal instability to be discussed 
in detail in Section VII, have been introduced for a gen- 
eral class of unbounded potentials including the PS as 
a particular case. This has been further elaborated and 
extended to include the GCM in Ref. 23. In both cases, 
the authors discuss directly the three-dimensional case, 
so that a direct comparison with the present work cannot 
be drawn at the present stage, but they also provide a 
detailed discussion of various approximations, within the 
general framework of density functional theory, that pro- 
vide a unified framework where even the present model 
could be included. 



V. EQUATION OF STATE 

Given an approximate solution of a fluid model there 
are several routes to the equation of state which, in gen- 
eral, give different results. The most common are three:^ 
the virial, the compressibility, and the energy route. The 
consistency of the outcome of these different routes can 
be regarded as an assessment on the soundness of the 
approximation. For some particular approximations it 
may also happen that the consistency of two of the three 
routes are automatically enforced (see Ref. 10 and ref- 
erences therein for a detailed discussion on this point). 
This is the case of the HPA, where the virial and energy 
routes coincide, as anticipated. Hence, the consistency 
with the compressibility route will provide a rough esti- 
mate of the regime of validity of the HPA within a phase 
diagram for the PSW potential. 



Let us briefly recall^ the methods to compute the com- 
pressibiUty factor, Z = jSp/ p^ associated with the three 
different routes. The virial route is defined by 

/>CXD 

Z = 1-/9/3/ dr ry (r) e-'''^^'' V (?-) , (5.1) 
Jo 

which, using standard manipulations,^ yields 

Z = l + p7,,[(l+7)y(l)-7Ay(A)]. (5.2) 

Thus the problem is reduced to the computation of the 
cavity function y{r) which, in the present context, follows 
fromEqs. (3.2) and (3.5). 

Next, we consider the compressibility route. 



z = i r^ 

P Jo S{0) ' 



(5.3) 



where the integral can be readily evaluated with the help 
of Eq. (3.9). 

Regarding the energy route, we start from the internal 
energy per particle 

1 /•! 

u = ^+er(l-7r)p / dry{r) 

-ea{l+la)pl dry{r). (5.4) 

In the above equation, the expressions given by Eqs. (2.1) 
and (3.1) have been used. In order to obtain Pp from u we 
exploit the following standard thermodynamic identity^" 



du\ _ fdpp\ 
dpjp ~ \dp)p 



thus leading to 



Z = \ 



P 



PI dp\^ I 

/o \dpj 13' 



du\ 



(5.5) 



(5.6) 



We have used the exact consistency between the virial 
and energy routes within the HPA as a test of the nu- 
merical calculations. 

Figure 1 depicts the results of the virial — and hence 
energy, as remarked — route under the condition e^/ca = 
2, which constitutes the borderline range of stability of 
the one-dimensional PSW model with A/cr < 1.^ Under 
this demanding condition, we have considered three re- 
duced temperatures from ksT/sa = 8 to ksT/sa ~ 12, 
whereas the width of the well has been fixed to the value 
A/cr = 1. We remark that pa is not limited in values 
from above due to the boundness of the potential. The 
clear downturn of all three curves for sufficiently large 
reduced density pa is a consequence of the existence of a 
maximum density pmax [see Eq. (BIO)], beyond which the 
HPA breaks down, as described at the end of Appendix 
B. In particular, the values of the maximum density for 
^r/^a = 2 and A/cr = 1 are Pmaxf = 4.94, 6.15, and 7.36 
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FIG. 1: Plot of the compressibility factor Z — Pp/p as 
computed from the virial(-energy) route to the pressure for 
Cr/ca ~ 2 and A/a = 1. Results for three different reduced 
temperatures ksT/ea = 8, 10, 12 are displayed. 



at kgT/ea — 8, 10, and 12, respectively, in agreement 
with Fig. 1. 

We compare in Fig. 2 the results from the virial (- 
energy) and the compressibility routes with MC 
simulations^^ for an intermediate value of the reduced 
temperature (fc^T/ca = 10) and other parameters as be- 
fore. Rather interestingly, the virial(-energy) route ap- 
pears to reproduce rather well the numerical simulation 
results up to the region where the artificial downward 
behavior shows up, whereas the compressibility route be- 
gins to deviate for densities pa > 3. 

As the temperature increases the HPA theory clearly 
remains a good approximation for a larger range of den- 
sities. We can naturally measure this by the requirement 
that virial(-energy) and compressibility routes are consis- 
tent within a few percent. This is indeed shown in Fig. 3, 
where we depict a transition line separating a "reliable" 
from an "unreliable" regime, as measured by the relative 
deviation of the two routes (here taken to be 5%), for four 
choices of the model parameters: (er/ea,A/cr) = (2,1), 
(2,0.5), (5, 1), and (5,0.5). The value A/cr = 0.5 is fre- 
quently used in the SW counterpart.^^ We observe that 
the region < p < p\im{T) where the HPA is reliable is 
hardly dependent on A/cr [compare curves (a)&(b) and 
(c)&(d) in Fig. 3]. On the other hand, at given val- 
ues of A/cr and ksT/ea, the range < p < pnmiT) 
decreases with increasing er/fa [compare curves (a)&(d) 
and (b)&(c) in Fig. 3], as expected. However, this effect 
is much less important if the increase of e^/ea takes place 
at fixed ksT/er (see inset of Fig. 3). It is interesting to 
note that, as illustrated by Fig. 2, the HPA virial route 
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FIG. 2; Comparison between the compressibility factor Z — 
Pp/a from the virial(-energy) and the compressibihty route 
to the pressure. The circles represent MC simulation results. 
Chosen parameters are er/ta = 2, A/a = 1, and ksT/ea = 
10. 



keeps being reliable up to a certain density higher than 

Aim- 

As said above, in Ref. 6 we introduced a LPA that was 
accurate for states where the penetrability effects were 
low or moderate. The LPA is complemented by the HPA 
presented in this paper. It is then interesting to compare 
the regions where each approximation can be considered 
reliable according to the same criterion as in Fig. 3. This 
is shown in Fig. 4 for the case e^/ea = 5 and A/cr = 1. 
The two transition lines split the plane into four regions: 
a region where only the LPA is reliable, a region where 
only the HPA is reliable, a region where both approxima- 
tions are reliable (and provide equivalent results), and a 
region where none of them is sufficiently good. The latter 
region shrinks as Sr/^a decreases (thanks to the HPA) 
or A/cr decreases (thanks to the LPA). 



VI. STRUCTURE 

As an additional test of the soundness of the HPA, we 
also study the RDF g{r) = h{r) + 1, which can easily be 
obtained from Eqs. (3.1) and (3.2), or equivalently from 
Eq. (3.6), once the auxiliary function w{r) has been de- 
termined. For a sufficiently high temperature (and hence 
high penetrability), the HPA is clearly well performing, 
as can be inferred from Fig. 5, when compared with stan- 
dard NVT MC results. Here we have considered the same 
parameters as in the preceding section (cr/ea = 2 and 
A/cr = 1) at a corresponding high-temperature value 




FIG. 3: Rough estimate of the region of reliability for the 
HPA, on the basis of the consistency between the virial(- 
energy) and compressibility routes, in the reduced density pa 
versus reduced temperature ksT/ea plane. Here the curves 
(a)-(d) correspond to {€r/€a,A/a) = (5,0.5), (5,1), (2,1), 
and (2, 0.5), respectively. The region below each curve repre- 
sents states where the relative deviation between the virial(- 
energy) route and the compressibility one is smaller than 5% 
and hence regarded as reliable. The inset shows the curves in 
the pa versus ksT/er plane for (cr/ta, A/a) = (2, 1) (dashed 
line) and (5,1) (solid line). 



ksT/ea = 10 and a density pa ~ 1.5 where overlap- 
pings are unavoidable. Under these conditions, there is 
no visible difference among the various approximations 
considered in Section IV. The excellent performance of 
the HPA observed in Fig. 5 agrees with the reliability cri- 
terion of Fig. 3 since the state /csT/ea = 10 and pa = 1.5 
is well below the curve (c) corresponding to e^/ca = 2 and 

A/cr=l. 

As we cool down, significant differences among var- 
ious approximations (HPA, nlHPA, mHPA, and MSA) 
begin to appear, as depicted in Fig. 6, where results cor- 
responding to temperatures ksT /ea = 5 (top panel) and 
ksT/ea = 3 (bottom panel) are reported within the same 
scale. The state ksT/ea = 5 and pa = 1.5 is still lying 
in the reliable region of Fig. 3, but close to the bound- 
ary line (c), while the state /cbT/eq = 3 and pa = 1.5 
is clearly outside that region. In the case ksT/ea = 5 
the HPA and its three variants are practically indistin- 
guishable, except in the region < r < a, which is very 
important to describe the correct thermodynamic behav- 
ior, where the best agreement with MC data corresponds 
to the nlHPA, followed by the HPA. The two approxi- 
mations that do not preserve the continuity of the cavity 
functions (mHPA and MSA) overestimate the jump at 
r = a. In the lower temperature case ksT/ea = 3 all 




FIG. 4: Regions of reliability for the LPA and the HPA in 
the reduced density pa versus reduced temperature ksT/ea 
plane. Here e^/ea^S and A/a = 1. The labels LPA, HPA, 
and LPA+HPA indicate the regions where only the LPA, only 
the HPA, or both approximations are reliable, respectively. 



the approximations overestimate the oscillations of the 
RDF. Interestingly, the HPA captures quite well the val- 
ues of g{r) near the origin. The worst overall behavior 
corresponds again to the MSA, which even predicts neg- 
ative values of g{r) for r/a < 1. 

Additional insights can be obtained by decreasing the 
range of interactions, in close analogy with what we con- 
sidered in previous work for the complementary LPA.^ 
This is reported in Fig. 7 for the case er/^a = 5 and 
A/a = 0.5 with {kBT/ea,pcr) = (5,1.5) (top panel) 
and {kBT/ca:PCT) = (2,0.8) (bottom panel). Again we 
stress that the same scale is used for both panels in 
order to emphasize the effect of lowering the tempera- 
ture. Clearly this is a more demanding situation. In 
fact, both states are above curve (a) of Fig. 3 and thus 
outside the corresponding reliability region. Therefore, 
clear deviations from MC results appear in all consid- 
ered approximations, especially in the lower temperature 
case ksT/ea — 2 (bottom panel). Yet, the HPA is still 
a reasonably good approximation that follows the main 
qualitative features of the correct g{r). In the higher 
temperature case ksT /ca = 5 the only noticeable limita- 
tions of the HPA practically take place near the origin, 
this deficiency being largely corrected by the mHPA. 

To conclude this section, it is worthwhile comparing 
the two complementary approaches HPA and LPA at a 
case where both are expected to be rehable, according to 
the diagram of Fig. 4. This is done in Fig. 8 for e^/ea = 5, 
A/cr = 1, ksT/ta — 2, and pa = 0.2. We observe that 
both approximations agree well each other and with MC 
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FIG. 5: Results for the RDF g(r) as a function of r/a with 
€r/ea = 2, A/a = 1, kaT/ea = 10, and pa = 1.5. Predic- 
tions from the HPA (solid line) are compared with MC results 
(circles). 



data, except in the region 1 < r/a < A/cr, where the 
HPA RDF presents an artificial curvature. 



VII. FISHER-WIDOM LINE AND 
FLUID-FLUID TRANSITION 



We now turn to an interesting point raised in previous 
work,® namely the question of whether the model can 
display a phase transition in spite of its one-dimensional 
character. 

The existence of general theorems — all essentially 
based on the original van Hove's result^^ — on the 
absence of phase transitions for a large class of one- 
dimensional models with short-range interactions is well 
established.^^ PSW and PS models, however, do not be- 
long to the class for which these general theorems hold. 
This is because boundness allows multiple partial (or 
even total) overlapping at some energy cost, thus render- 
ing invalid the arguments used in the aforementioned 
theorems. 

On the other hand, none of these theorems pro- 
vide a general guideline to understand whether a one- 
dimensional model may or may not display a non-trivial 
phase transition, and one has then to rely on the speci- 
ficity of each model. 

As discussed in our previous work,® it is instructive 
to first address the simpler question of the location of 
the Fisher-Widom (FW) line. This is a line separat- 
ing two different regimes for the large-distance behav- 
ior of the RDF g{r) in the presence of competing repul- 
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FIG. 6: Comparison of different approximations in the results 
for the RDF g{r) vs r/a. Chosen parameters are er/ta = 2, 
A/a = 1, and pa = 1.5. MC results (circles) are compared 
with the HPA (solid line), the nlHPA (dashed line), the mHPA 
(short dashed line), and the MSA (long dashed line). The 
two different panels refer to different reduced temperatures: 
ksT/ea = 5 (top panel) and ksT/e^ = 3 (bottom panel). 
Note that both panels are drawn within the same scale. 



FIG. 7: An additional comparison of different approximations 
in the results for the RDF g(r) vs r/a. Here fixed parameters 
are er/fa = 5 and A/cr = 0.5. MC results (circles) are com- 
pared with the HPA (solid line), the nlHPA (dashed line), the 
mHPA (short dashed line), and the MSA (long dashed line). 
The top panel refers to the state ksT/ta = 5 and pa = 1.5, 
whereas the bottom panel refers to ksT/ta = 2 and pa = 0.8. 
In the bottom panel the MSA is not depicted. Again both 
panels are on the same scale. 
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FIG. 8: Comparison of different approximations in the results 
for the RDF g{r) vs r/a. Chosen parameters are tr/ea = 5, 
A/a = 1, ksT/ea = 2, and pa = 0.2. MC results (circles) 
are compared with the HPA (solid line) and the LPA (dashed 
line). 

sive/attractivc interactions.^^ The rationale behind the 
FW line is that on approaching the critical points where 
attractions become more and more effective, the behav- 
ior of correlation functions must switch from oscillatory 
(characteristic of repulsive interactions) to exponential 
with a well defined correlation length ^. In the previ- 
ous work, we analyzed the location of this line for PSW 
within the LPA. Here we extend this analysis to the HPA 
regime and discuss the compatibility of the two results. 

Let us first briefly recall the main points of the analysis, 
referring to Ref. 6 for details. From Eq. (3.1) we note 
that the asymptotic behavior of g{r) is the same as that 
of y{r). In view of Eq. (3.2), this is hence related to 
w{r), whose asymptotic behavior is governed by the pair 
of conjugate poles of w{k) with an imaginary part closest 
to the origin. If the real part of the pair is zero, the decay 
is monotonic and oscillatory otherwise. 

According to Eq. (3.3), the poles of w{k) are given by 



PlrfsW (k) 



1. 



(7.1) 



Let k = Hx be the imaginary pole and k = ±i(x' ± iy) 
be the pole with the imaginary part closest to the origin. 
The FW line is determined by the condition x = x'. 
This gives, at a given temperature, three equations in 
the three unknowns x, y, and p.^ More specifically, after 
some algebra, one gets 

y . . 

— sinhx sinh(Ax) [cosy — cos(Ay)] = sinha;cosh(A.T) 

X 

X sin(Ay)(cos J/ — 1) — sinh(Aa;) cosh a; 
X siny [cos(Ay) — 1] , (7.2) 



sinh(Aa::) cos(Ay) — 1 



n -1 



sinh ; 



cos y — 1 



-1 



1 



P = 



27r 7sinh(Ax') — (1 + 7) sinhx 



(7.3) 



(7.4) 



The inverse of the parameter x represents the correla- 
tion length ^ = 1/x. From a practical point of view it is 
more convenient to use x rather than ksT /ca as a free 
parameter to construct the FW line. In that way, Eq. 

(7.2) becomes a transcendental equation that gives y as 
a function of x\ once y{x) is known, the solution to Eq. 

(7.3) gives ksT/ea as a function of x; finally, insertion 
of y(x) and kBT{x)/ta into Eq. (7.4) provides p{x). The 
corresponding values of the pressure are obtained from 
either Eq. (5.2) (virial-energy route) or Eq. (5.3) (com- 
pressibility route). 

We observe that T decreases as x decreases, until a 
critical value Tc is found in the limit x — )• 0. In that 
limit, Eqs. (7.2)-(7.4) simplify to 



Aye [cos yc - cos(Ayc)] = 



sin( Aye) (cos yc - 1) - A sin yc 
X [cos(Ayc) - 1] , (7.5) 



cos (Aye) - 1 



cos yc 



1 



1 



27^ 7cA - 1 



(7.6) 



(7.7) 



At the critical point T = Tc and p = pc^ one has x = or, 
equivalently, ^ ^- 00. Therefore, at this point w{r) does 
not decay for long distances and in Fourier space one has 
w{k) ^ fc~^ and S{k) ^ fc~^ for short wave numbers. 
The condition S{k) ^ fc^^ is also satisfied for T < Tc 
if pa = [27r(7A/o' — 1)] , in agreement with Eq. (3.9). 
This defines in the p-T plane a spinodal line or locus of 
points of infinite isothermal compressibility (within the 
compressibility route). The spinodal line cannot be ex- 
tended to temperatures larger than the critical value Tc 
because [27r(7A/o'— 1)] > Pmaxf if T > Tc, where 
Pmax is the maximum density beyond which the HPA 
is unphysical at a given temperature (see Appendix B). 
We further note that the spinodal line only has a lower 
density (or vapor-like) branch, thus hampering the inter- 
pretation of {Tc, Pc) as a conventional critical point. 

The above features are already suggestive of consider- 
ing the HPA spinodal line as an artifact of the theory 
when used in a region of parameter space where the ap- 
proximation is invalid. Additional support to this view 
stems from the fact that the HPA keeps predicting a 
spinodal line and a critical point even in the SW case 
[trlf-a — !• 00, 7,- — !• 1), a clearly incorrect feature. As we 
shall discuss further below, specialized numerical simula- 
tions coupled with a recent analytical study, ^^ strongly 
support the absence of any phase transition in the present 
one-dimensional PSW model. 
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plane with A/cr = 1 and er/fa = 10 (top panel) and e^/ea = 5 
(bottom panel) . The long dashed curves are the spinodal lines 
predicted by the HPA. The FW and spinodal liness meet at 
the critical point (denoted by a circle). 
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In Fig. 9 we report the comparison between the FW 
lines, as predicted by the LPA and the HPA, in the pa vs 
ksT/ca plane for A/tr = 1 and two different energy ratios 
^r/^a = 10 and er/ea = 5. The spinodal line predicted 
by the HPA is also included. As said above, the FW and 
spinodal lines meet at the critical point. While at high 
temperatures (above fcsT/ea = 3) there is a remarkable 
agreement between the two approximations, deviations 
occur at lower temperatures. 

As in the original work by Fisher and Widom,^^ we 
also report the FW line in the pu/ta vs ksT/ea plane 
(see Fig. 10), where again we compare the lines as de- 
rived from the HPA and LPA schemes. For the energy 
ratio Er/^a = 10 (top panel) we see that the HPA and 
the LPA give qualitative similar forms of the FW line 
with a significant deviation at low temperatures, where 
again the HPA FW line is interrupted at T = T^- Note 
that, while all three standard routes give practically iden- 
tical results within the LPA, the virial(-energy) route 
in the HPA differs from the compressibility one at low 
temperatures (more than 5% for ksT/ea ^ 1.2). For 
consistency, the HPA spinodal line is obtained via the 
compressibility route only. Similar features occur for the 
second lower value of the energy ratio, namely er/ea = 5 
(bottom panel) . Again, the LPA and HPA lines are qual- 
itatively similar, the three routes in the LPA provide in- 
distinguishable results, and the virial and energy routes 
in the HPA deviate more than 5% for ksT/ea < 1.2. The 
main distinctive feature in this case er/ea = 5 is a marked 
upswing of the tail of the FW line, absent the previous 
case er/ea = 10. This means that, on increasing pene- 
trability — that is, on decreasing er/ea — the transition 
from oscillatory (above the line) to monotonic (below the 
line) behaviors occurs at a higher pressure and a higher 
density for a fixed temperature ksT/ea- 

Despite the important differences in the steps followed 
to derive the LPA and the HPA, it is noteworthy that 
they agree in the qualitative shape of the FW lines (even 
though the HPA predicts a spurious spinodal line). It is 
then reasonable to expect that the true FW line should 
interpolate the LPA line at low temperatures with the 
HPA line at high temperatures. 

Finally, we now tackle the issue of the existence of a 
phase transition for the PSW one-dimensional model. In 
view of the HPA results on the seemingly existence of 
a spinodal line (and hence of a critical point), we here 
consider the fluid-fluid transition. As we shall see, our 
numerics is compatible with the absence of such a transi- 
tion, thus supporting the view that the above findings of 
a spinodal line is indeed a consequence of the application 
of the HPA to a regime where the theory is not valid. 

As the FW line always anticipates the critical point, as 
remarked, we can then look for the existence of a fluid- 
fluid coexistence line in the region predicted by the in- 
terpolation of the LPA and HPA Fisher-Widom lines. 
We have carried out extensive simulations of the PSW 
fluid using Gibbs ensemble Monte Carlo techniques and 
employing all standard improvements suggested in the 
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FIG. 10: Same as in Fig. 9, but in the pcr/ea vs kBT/ea 
plane. Note that while in the LPA the three routes to the 
pressure are not distinguishable one from the other on the 
graph scale, for the HPA the difference between the virial and 
the compressibility route is noticeable at low temperatures. 



literature. ^'^"^^ In order to validate our code, we tested 
it against the case of the one-dimensional SW potential, 
where exact analytical predictions for all thermodynamic 
quantities are available. 

We have used up to 1000 particles and carefully 
scanned the temperature range 0.1 < ksT/ea < 2.0 and 
the density range 0.1 < per < 6, as suggested by the FW 
line (see Fig. 9). We have also considered different values 
of er/^a and A/a for cases giving a significant overlap- 
ping probability. In all the cases we have not found any 
signature of a fluid-fluid phase separation. 

Although the absence of a critical transition is always 
much more difficult to assess as opposed to its presence, 
the first scenario is consistent with more than one indi- 
cation. The first indication stems from a lattice model 
counterpart of the one-dimensional PSW model. This is 
discussed in Appendix C, where the lattice version of the 
PSW model is constructed following standard manipula- 
tions with the result that no phase transitions are present 
for finite occupancy. An additional evidence supporting 
the absence of any fluid-fluid or freezing transition stems 
from the very recent exact analytical work alluded to 
earlier, ^^ that, using the methodology presented in Ref. 
18, concludes that no phase transitions are present for 
the PSW and PS models in one dimension. 

In our simulations we have also investigated values of 
Cr/ea and A/a which violate Ruelle's stability criterion 
(see Appendix A) and thus the one-dimensional PSW 
model is not necessarily stable in the thermodynamic 
limit. Here the phenomenology turns out to be much 
more interesting. For sufficiently low temperatures and 
sufficiently high densities we observe the formation of a 
"blob" of many-particle clusters (each made of a large 
number of overlapping particles) having a well defined 
and regular distribution on the axis, and occupying only 
a portion of the system length. In the "blob" phases 
the energy per particle grows with the number of par- 
ticles, thus revealing the absence of a thermodynamic 
limit. The transition from a "normal" phase to a "blob" 
phase if e^/ea < 2(£ -I- 1), where £ is the integer part of 
A/ a, is illustrated in Fig. 11. This figure shows the RDF 
at per = 6 (top panel) and pa ^ 8 (bottom panel) for 
e^/ea = 2, A/cr = 1.5, and ksT/ca = 10. At the lower 
density the structure of the PSW fluid is qualitatively 
not much different from that expected if Er/ea > 2{£+ 1) 
(compare, for instance, with Fig. 6). However, the struc- 
ture at the higher density is reminiscent of that of a 
solid, except that the distribution of particles does not 
span the whole length (here L — N/p = 62.5cr). In- 
stead, the particles distribute into a few clusters regu- 
larly spaced with distance 1.17a, so that the particles of 
a given cluster interact attractively with all the particles 
of the nearest- and next-nearest-neighbor clusters. Note 
that the number of particles within any given cluster is 
not necessarily identical. It is also worth noticing that, in 
spite of the huge difference in the vertical scales of both 
panels in Fig. 11, they are consistent with the condition 
J^^/'drh{r)^-l/2p. 
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FIG. 11: Plot of the RDF g{r) obtained from MC simulations 
for er/ta = 2, A/a — 1.5, ksT/ea = 10, and pa — 6 (top 
panel) and pa — 8 (bottom panel). 



The decay of the peaks oi g{r) is mainly due to the lack 
of translational invariance, i.e., the first and last clusters 
have only one nearest-neighbor cluster, the first, second, 
next-to-last, and last clusters have only one next- nearest- 
neighbor cluster, and so on. 

We stress that the above phenomenon is specific of 
bounded potentials, such as PSW, and has no counter- 
part in the hard-core domain. It is then plausible to 



expect their appearance even in the corresponding three- 
dimensional versions of these models where freezing tran- 
sition (and phase separation for SW) are present but 
could both be hampered by the presence of this clustering 
phenomenon in the region of parameter space where Ru- 
cUe's stability criterion is violated. This would extend the 
interesting phenomenology already established for the PS 
case^^. Work along these lines on the three-dimensional 
case is underway and will be reported elsewhere. 



VIII. CONCLUSIONS 

In this paper we have completed the study initiated 
in previous work^'^ on the penetrable-square-well model 
in one-dimension. This is a model combining the three 
main ingredients present in many colloidal-polymer solu- 
tions, namely repulsions, attractions, and penetrability. 
While the first two are ubiquitous even in simple fluids, 
the latter is a peculiarity of complex fluids where there 
exist many examples of colloid-polymer systems which 
are penetrable (with some energy cost) to some extent, 
and they involve both steric repulsions and short-range 
attractions. This model then captures all these crucial 
features at the simplest level of description within an im- 
plicit solvent description. 

The main new point of this paper was to present an ad- 
ditional and complementary approximation, denoted as 
the HPA, valid in regimes complementary to those valid 
for the low-penetrability scheme LPA discussed in Ref. 6. 
While the idea behind the LPA was to modify the exact 
relations valid for the one-dimensional SW fluid — and 
in fact for any fluid with a hard core and short-range 
attraction — to allow penetrability within some reason- 
able approximation, the driving force behind the HPA is 
the fact that for bounded potential the Mayer function 
can be made arbitrarily small by considering sufficiently 
high temperatures. As a consequence, only the linear 
chain diagrams need to be retained at each order in the 
cluster expansion. As it turns out,^° this is tantamount 
to considering the celebrated Dcbye-Hiickel theory for 
charged fluids, and we have here considered the sound- 
ness of this approximations at various regimes as com- 
pared to specialized MC simulations. The latter were 
also compared with other approximations which parallel 
the entire hierarchy of approximations in the framework 
of charged fluids, ranging from the most sophisticated 
non-linear Debye-Hiickel to the simplest mean-spherical 
approximations. We have assessed the regime of reliabil- 
ity of these approximations both for thermodynamic and 
correlation functions by comparison with MC simulations 
and by internal consistency between different routes to 
thermodynamics. 

Next we have also discussed the location of the FW 
line, separating oscillatory from monotonic behavior in 
the correlation function, within the HPA and compared 
with that obtained from the LPA introduced in previous 
work.^ In agreement with previous findings, we find that 
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penetrability enhances the region where eorrelation func- 
tions have a monotonic regime. The FW derived from 
the HPA and LPA schemes are found to be in quahtative 
agreement thus making it possible the drawing of a line 
interpolating the high- and low-temperature regimes. 

As a final point, we investigated the possibility of a 
fluid-fluid transition. This possibility arises because the 
boundness of the potential renders the van Hove theorem 
on the absence of phase transition for one-dimensional 
model with short-range interactions non-applicable. In 
fact, the HPA is seen to predict a critical point where 
the FW line meets a spinodal line. However, this predic- 
tion takes place in a region of densities and temperatures 
where the HPA is not reliable. A careful investigation 
using both NVT and Gibbs ensemble Monte Carlo tech- 
niques akin to those exploited in the investigations of 
the analog problem for the three-dimensional SW model 
yields negative results. These findings are also supported 
by analytical arguments based on the lattice gas coun- 
terpart where the absence of transition can be motivated 
by the absence of an infinite occupancy of each site, as 
well as by an exact analytical proof^^ of a no-go theo- 
rem proving the absence of any phase transition in this 
model, which is in agreement with, and beautifully com- 
plements, our work. 

In our quest for a possible thermodynamic transition in 
the one-dimensional PSW model we have explored values 
of the energy ratio e,./ea and well width A/a for which 
the stability of the system in the thermodynamic limit is 
not guaranteed by Ruelle's criterion. ^^'^^ We have found 
that, as the temperature decreases and/or the density 
increases, a transition from a normal fiuid phase to a 
peculiar solid-like phase takes place. The latter phase 
is characterized by the formation of clusters of overlap- 
ping particles occupying a small fraction of the available 
space and with non-extensive properties. This clustering 
transition preempts both the fiuid-fiuid and fiuid-solid 
transitions. 

In view of the results presented here, it would be very 
interesting to discuss the phase diagram of the corre- 
sponding three-dimensional PSW model. The phase dia- 
gram of the SW model {er/ea ~^ oo) is indeed well estab- 
lished and includes both a fluid-solid transition — present 
even in the HS counterpart — and a fluid-fluid transition 
line. The latter is present for any value of ea and A/ a but 
is stable against freezing only for A/a > 0.25, the depth 
of the well being irrelevant. ^^ The results presented here 
strongly suggest the importance of the additional param- 
eter ir/^a- A first interesting issue would be the Ruelle 
instability in three-dimensions. A straightforward exten- 
sion of the arguments presented in Appendix A predicts a 
guaranteed stability for e^/ea > 12 (if A/a < 1), but the 
actual onset of the instability cannot be assessed through 
these arguments. One could expect that for sufficiently 
high penetrability (i.e., er/ca < 12 and high density) 
a phenomenon akin to the "clustering" transition found 
here could be present. In the case Cr/ca > 12, where the 
clustering transition is not expected, the interesting point 



is to assess the influence of the ratio Cr/ea on the location 
of the fluid- fluid critical point and coexistence line. All of 
this opens the possibility of a rich and interesting phase 
diagram which would complement that already present in 
a general class of bounded potentials with no attractive 
tails. ^^'^'^ We note that the high-penetrability regime is 
indeed the realm of the HPA presented here, which can 
be obviously extended to three-dimensions. 

Work on the three-dimensional PSW model including 
the above points and other aspects is underway and will 
be reported elsewhere. 
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Appendix A: Ruelle's stability criterion 

Let us consider the ID PSW model characterized by 
er/ea and A/a. Let us call i the integer part of A/a, 
i.e., £ < A/ a < £ + 1. According to Ruelle's criterion, a 
sufficient condition of thermodynamic stability is^^'^^ 

N-l N 

Un{xi,...,xn)=Y1 Y. <t>{\x^~XJ\)>-NB (Al) 

Z— 1 J— 2 + 1 

for all configurations {xi}, where i? is a fixed bound. 

Given the number of particles N, we want to obtain the 
configuration with the minimum potential energy Un. 
Without loss of generality we can see any given configu- 
ration as a set of M clusters (1 < M < N), each cluster 
i being made of Sj overlapping particles (i.e., any pair 
of particles of a given cluster are separated a distance 
smaller than a) . In Rcf. 5 we proved that for a fixed value 
of M the minimum energy corresponds to Si = s = N/M, 
all the particles of each cluster being located at the same 
point and the centers of two adjacent clusters being sep- 
arated a distance a. Therefore, we can restrict ourselves 
to this class of ordered configurations and use s as the 
variational variable. 

The repulsive contribution to the potential energy is 



U^^is) = M 



s{s~l) 



(A2) 



To compute the attractive contribution we need to take 
into account that all the particles of a given cluster in- 
teract attractively with the particles of the £ +1 nearest 
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clusters. The total number of pairs of interacting clusters We then see that the value that minimizes Un{s) is 
are {e + l)[M -{£+!)]+ £+{£-!) + {£ -2) + ■■■ + 1 = 
{£ + 1) (Af - 1 - £/2). Therefore, 



UUs) = -{e+l)[M-^-±^].s^e„ 



The total potential energy Un ~ UJ^ + W^ is 



(A3) 



Un{s) = N^er 



S^ 



N 
1+2 



1 



2i£+l)ea 



(A5) 



{£+l)[N-^-^s]sea. (A4) ofC/wis 
I 



This value is only meaningful if er/2{£ + l)ea < 1- Oth- 
erwise, s* = 1. In summary, the absolute minimum value 



U*^ = Un{s.) = 



-f^e„ 



-(^ + l)(iV- 






i2 



2(«+l)c 



er<2{e+l)ea, 
er>2{£+l)ea. 



(A6) 



Therefore, if e,, > 2{£ + l)ea the potential energy is 
bounded from below by ~NB with B = {£+l)ea and thus 
the system is stable in the thermodynamic limit. On the 
other hand, if e^ < 2{£+ l)ea there exist configurations 
that violates Ruelle's criterion and so the thermodynamic 
stability of the system is not guaranteed. 



Equation (3.4) can be rewritten as 



/sw(fc) = T [(1 + 7) (e*' - e-'') - 7 (e''^ - e~^''^)] 



Therefore, 



(B4) 



Appendix B: Density expansion of w{r) within the 
HPA 



Starting from Eq. (3.3) and for P7r|/sw(^")l < 1; the 
Fourier transform w[k) can be expanded in power series 
as 



w{k) = ^(p7.)""V^W(fc) ■ 
Upon inverse Fourier transform one then has 

oo 
w{r) ^ ^ (P7r)"~^ Wn {r) , 

where 



)i=2 



, Zn 



(Bl) 



(B2) 



(B3) 



.„ n 771 n — -rn / \ / \ / 

m=Op=0 q=0 \ / \^/ \ y 
X(-1)™+P+9(1+^)™^"-™ 
y,^ik[2p-m+{2q-n+m)X] (B5) 

The origin (fc = 0) is a regular point of fsw{k) and 

hence of /sw(^) [^^^ ^'^^ ^^ each separate term in Eq. 
(B5)], so we can choose to save the point fc = in Eq. 
(B3) either from above or from below. Here we do it from 
above with the result 



hm ^- I *^^ 

e^o+ Jc 27r {k + icy 



i-ry 



(n-1) 



-e(-r), n>0, 



(B6) 

where the path C in the complex k plane goes from k = 
— cx) to fc = -|-cx) and closes itself on the upper plane if 
r > and in the lower one if r < 0. In Eq. (B3) we then 
find 



IL in It — III / \ / \ / \ 



\—r — 2p + m — {2q — n + m) A] 



n-l 



-6 [-r -2p + m- {2q-n + m) A] 



(B7) 



It is interesting to note that io„(r) = if r > An. Thus Eq. (B2) can be rewritten as 



n— niax{2,[r/A] } 
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> 

o 

o 

Q. 

?^ 
D 



?>- 



2.5 



2.0 



/Tr Pmax^ (A/a=0.5) 

-Tr Pconv^ (A/a=0.5) 
-Tr Pconv^ i^^=^ -0) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

yA/a 

FIG. 12: Plot of the radius of convergence 7rPconvO" and of 
the maximum density 7rPmaxO" versus ^A/a for A/a = 0.5 
and A/cr = 1. 



Therefore, w(r) is still well defined beyond the radius of 
convergence, i.e., for p > Pconv On the other hand, anal- 
ogously to the PS case,20.29 the HPA for the PSW fluid 
becomes unphysical, at a given temperature, for densities 
larger than a certain value Pmax given by the condition 



yfrPmax) 



f<. 



SW 



(BIO) 



where f^^^ is the absolute maximum value of /sw(fc)- 



Since f. 



/■max 
/SW 
max <-- 

SW -i 



SW , it is obvious that p^ax > Pconv 

max 

For sufhciently large values of 7A (actually, for temper- 
atures below the critical value Tc defined in Section VII) 

one has f^"" = /sw(0), and so jrPnia,^ = (7A - 1)"V2- 
In that case, the line of maximum density becomes a 
spinodal line, as discussed in Section VII. Figure 12 also 
includes a plot of 7rPmax as a function of 7A for the 
same two values of A/cr. Note the kink of the curve 
TrPmax for A/a = 1 at the critical point (-fc^/cr ~ 2.11, 
7rPc — 0.45), so that 7rPmax = (7A — l)^"^/2 (spinodal 
line) if 7 A/cr > 2.11. 



Appendix C: The PSW lattice gas 



where [r/X] is the integer part of r/X. 

The radius of convergence of the series (B8) depends on 
temperature and can be obtained by the same arguments 
as in the PS case.^ From the denominator of Eq. (3.3), it 
follows that the series converges provided that p < Pconv, 
where 



\^r Pconv ) 



fsw 



(B9) 



Here 



SW 



denotes the absolute maximum value of 



/sw(fc) • From Eq. (3.4) if 7 < (A^ - 1)-^ that max- 



imum corresponds to k = 0, i.e.. 



SW 



-/sw(0), 



(1 - 7A)"V2- On the other hand. 



/sw 



takes 



and so jrPcom 

if 7 > (A"^ — 1)~^ the maximum value 

place at fc 7^ and so 7rPconv must be obtained nu- 
merically. For sufficiently large values of 7A one has 

/sw = /sw(0), so that 7rPconv = (7^ - 1)"V2 and 

max 

this coincides with the maximum physical density (see 
below). Figure 12 shows ^rPconv as a function of 7 A for 
two values of A/cr. In the PS limit (7A — > 0) one has 
JrPconv = 2 ■ ^'^ ^^'^ strength of the attractive part of the 
potential (measured by the product 7A) increases, the 
radius of convergence first grows, reaches a maximum, 
and then decays. 



Except when 7A is so large that 



./s 



w 



= /sw(0), 



the maximum 



/sw 



corresponds to a negative value 



of /sw(^) and so the singularity responsible for the ra- 
dius of convergence is located on the negative real axis. 



Consider the PSW model in one dimension. The grand 
partition function is 



1 /e'^'' 



N=0 



s(p,i,r) = ^— — z^(i,r), (ci) 



N 



where A^ is the thermal dc Broglie's wavelength and 



Zn {L,T) = / dri . . . drjv exp 



N 






'(r« 



(C2) 

We now discretize the length L as a sum of A^c ^ 1 cells 
of size a = L/Nc with occupancy tt-q,, a ~ l,...,Nc. 
The value of a is chosen in the interval a < a < a + 
A, with A < cr, so that two particles in the same cell 
are assumed to interact repulsively and two particles in 
adjacent cells are assumed to interact attractively. The 
integral can then approximated as J, dr « aNc and the 
configurational partition function Zjy as 



Zn{L,T) « a^ ^ (SAT^^^n^exp 

xA^!. 



1 ^' 



al3=l 



(C3) 



In Eq. (C3) the 5 function accounts for the constraint 
A^ = ^^ Ua and the factor A^! of the indistinguishability 
of the particles. Also we have introduced Xq/3 accounting 
for multi-interactions among cells and 0^/3 which is equal 
to £,, ii a — j3 and equal to — e^ for nearest-neighboring 
(nn) cells. Noting that each particle in an Q!-cell can 
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either interact with n^ — 1 other particles within the same 
cell or with np particles within a nn /3-cell, we see that 



Substituting into Eq. (CI) we then find 



Afe 



9 2^ Xa/30c</3 



Q,/3 = l 



Ea ^ naUfj + -tr ^ Tl^ («„ - 1) 



(«,^> 



(C4) 



H(M,L.r, . Y. il 



N=0 



=./3m 



Af 






4 Ex. 



^0Q,3 



Q,^ 



(C5) 



Because of the S, the two sums can be inverted and the sum over A^ can be explicitly carried out thus obtaining 



S (/i, L,T) = ^ cxp 

{n„} 






-er ^ ?i„ (na - 1) - /i ^ Wq 



(C6) 



with /i = M + (1//3) ln(a/AT). 

Assume a finite length L (and hence a finite number of cells Nc) with periodic boundary conditions. The above 
partition function can then be solved by standard transfer matrix techniques 



Na 



{ni ....TiiVc } '^ — 1 



where we have introduced the matrix 

iT-an^ — exp 



/? ( -eanaUp + -J {Ua {Ua ~ I) + Tip {np - 1)) - - (^a + TLp) 



(C7) 



(C8) 



r 



If Nc is finite, one then has in the thermodynamic limit 
hm ilogS(/i,i,r) = ilogAo, (C9) 



where Aq is the largest eigenvalue of the matrix. Clearly 
this is analytic and no phase transitions are possible for 
finite occupancy N^- 
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